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Abstract 

The growing complexity of real-world problems has motivated computer scientists to search for efficient problem-solving 
methods. Metaheuristics based on evolutionary computation and swarm intelligence are outstanding examples of nature- 
inspired solution techniques. Inspired by the social spiders, we propose a novel Social Spider Algorithm to solve global 
optimization problems. This algorithm is mainly based on the foraging strategy of social spiders, utilizing the vibrations 
on the spider web to determine the positions of preys. Different from the previously proposed swarm intelligence 
algorithms, we introduce a new social animal foraging strategy model to solve optimization problems. In addition, we 
perform preliminary parameter sensitivity analysis for our proposed algorithm, developing guidelines for choosing the 
parameter values. The Social Spider Algorithm is evaluated by a series of widely-used benchmark functions, and our 
proposed algorithm has superior performance compared with other state-of-the-art metaheuristics. 
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1. Introduction 

W ITH the fast growing size and complexity of modern 
optimization problems, evolutionary computing is 
becoming increasingly attractive as an efficient tool for op¬ 
timization. Depending on the nature of phenomenon simu¬ 
lated, evolutionary computing algorithms can be classified 
into two important groups: evolutionary algorithms (EAs) 
and swarm intelligence based algorithms. EAs, which mainly 
draw inspiration from nature, have been shown to be very 
successful for optimization among all the methods devised 
by the evolutionary computation community. Currently 
several types of EAs have been widely employed to solve 
real world combinatorial or global optimization problems, 
including Genetic Algorithm (GA), Genetic programming 
(GP), Evolutionary Strategy (ES) and Differential Evolu¬ 
tion (DE). These algorithms demonstrate satisfactory per¬ 
formance compared with conventional optimization tech¬ 
niques, especially when applied to solve non-convex opti¬ 
mization problems m- 

In the past two decades, swarm intelligence, a new kind 
of evolutionary computing technique, has attracted much 
research interest [5]. The term swarm is employed in a 
general manner to refer to any collection of interactive 
agents. Swarm intelligence is mainly concerned with the 
methodology to model the behavior of social animals and 
insects for problem solving. Researchers devised optimiza¬ 
tion algorithms by mimicking the behavior of ants, bees, 
bacteria, fireflies and other organisms. The impetus of cre¬ 
ating such algorithms was provided by the growing needs 
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to solve optimization problems that were very difficult or 
even considered intractable. 

Among the commonly seen animals, spiders have been 
a major research subject in bionic engineering for many 
years. However, most research related to spiders focused 
on the imitation of its walking pattern to design robots, 
e.g. [4]. A possible reason for this is that a majority of 
the spiders observed are solitary [5], which means that 
they spend most of their lives without interacting with 
others of their species. However, among the 35 000 spider 
species observed and described by scientists, some species 
are social. These spiders, e.g. Mallos gregalis and Oecobius 
civitas, live in groups and interact with others in the same 
group. Based on these social spiders, this paper formulates 
a new global optimization method to solve optimization 
problems. 

Spiders are air-breathing arthropods. They have eight 
legs and chelicerae with fangs. Spiders have been found 
worldwide and are one of the most diverged species among 
all groups of organisms. They use a wide range of strate¬ 
gies for foraging, and most of them detect prey by sensing 
vibrations. Spiders have long been known to be very sen¬ 
sitive to vibratory stimulation, as vibrations on their webs 
notify them of the capture of prey. If the vibrations are in 
a defined range of frequency, spiders attack the vibration 
source. The social spiders can also distinguish vibrations 
generated by the prey with ones generated by other spi¬ 
ders [6] . The social spiders passively receive the vibrations 
generated by other spiders on the same web to have a clear 
view of the web. This is one of the unique characteristics 
which distinguishes the social spiders from other organisms 
as the latter usually exchange information actively, which 
reduces the information loss to some degree but increases 
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the energy used per communication action [7]. 

The group living phenomenon has been studied inten¬ 
sively in animal behavior ecology. One of the reasons that 
animals gather and live together is to increase the possi¬ 
bility of successful foraging and reduce the energy cost in 
this process [H]. In order to facilitate the analysis of so¬ 
cial foraging behaviour, researchers proposed two foraging 
models: information sharing (IS) model [5] and producer- 
scrounger (PS) model [TU]. The individuals under the IS 
model perform searching and seek for opportunity to join 
other individuals simultaneously. In the PS model, the 
individuals are divided into leaders and followers. Since 
there is no leader in social spiders m, the IS model is 
more suitable to formulate the foraging behavior of social 
spiders, and we use this model to control the searching 
pattern of our proposed algorithm. 

In this paper, inspired by the social behavior of the so¬ 
cial spiders, especially their foraging behavior, we propose 
a new metaheuristic for global optimization: the Social 
Spider Algorithm (SSA). The foraging behavior of the so¬ 
cial spider can be described as the cooperative movement 
of the spiders towards the food source position. The spi¬ 
ders receive and analyze the vibrations propagated on the 
web to determine the potential direction of a food source 
[12] . We utilize this natural behavior to perform optimiza¬ 
tion over the search space in SSA. 

The contribution of this paper is threefold: 

• We propose a new nature-inspired swarm intelligence 
algorithm based on social spiders. This population- 
based general-purpose metaheuristic demonstrates out¬ 
standing performance in the global optimization bench¬ 
mark tests. 

• We introduce a new social animal foraging model 
into metaheuristic design. This is the very first at¬ 
tempt of employing the IS model to solve optimiza¬ 
tion problems. We also incorporate the information 
loss schemes in the algorithm, which is a unique de¬ 
sign of our proposed algorithm. 

• We perform a series of experiments to investigate the 
impact of different parameters and searching schemes 
on the performance of the algorithm. The result of 
these experiments may serve as important inputs for 
further research. 

The rest of this paper is organized as follows. We will 
first present some related work on swarm intelligence and 
bio-inspired metaheuristics in Section]^ Then we will for¬ 
mulate and elaborate on SSA by idealizing and imitating 
the foraging behavior of social spiders in Sectionj^ Section 
1^ introduces the benchmark functions we use for testing 
the performance of SSA, with the experimental settings. 
Section presents the simulation results of SSA on the 
benchmark functions and the comparison with other pop¬ 
ular metaheuristics. Finally we will conclude this paper in 
Section]^ and propose some future work. 


2. Background 

Swarm intelligence algorithms mimic the methods in 
nature to drive a search for the optimal solution. At the 
very beginning there are two major methods for this kind 
of algorithms: ant colony optimization (AGO) [13] and 
particle swarm optimization (PSO) |I4j . 

AGO is inspired by the foraging behavior of ants, whose 
goal is to find a shortest path from their colony to food 
sources. In this metaheuristic, feasible solutions of the 
optimization problem to be solved are represented by the 
paths between the colony and food sources. The ants com¬ 
municate with and influence others using pheromone, a 
volatile chemical substance. When an ant finds a food 
source, it deposits certain amount of pheromone along the 
path and the amount is positively correlated with the qual¬ 
ity of the food source. The pheromone laid down biases the 
path selection of other ants, providing positive feedback. 
Using the scheme of positive feedback, the algorithm leads 
the ants to find the shortest path to a best food source 

m- 

PSO is motivated by the movement of organisms as a 
group, as in a flock of birds or a school of fishes. The 
group is represented by a swarm of particles and PSO uses 
their positions in the search space to represent the feasi¬ 
ble solutions of the optimization problem. PSO manipu¬ 
lates the movement of these particles to perform optimiza¬ 
tion, utilizing the information of individual experience and 
socio-cognitive tendency. These two kinds of information 
correspond to cognitive learning and social learning, re¬ 
spectively, and lead the population to find a best way to 
perform optimization |14] . 

The above two metaheuristics have been applied to 
solve a vast range of different problems, e.g. [I5|[T6]. Mo¬ 
tivated by such success, swarm intelligence algorithm de¬ 
sign has attracted many researchers and several new algo¬ 
rithms were devised. The most widely studied organism in 
swarm intelligence is the bee |3|. Abbass proposed a Mar¬ 
riage in honey Bees Optimization (MBO) in [T^ and this 
algorithm was applied to solve propositional satisfiability 
problems (3-SAT problems). In MBO, the mating flight of 
the queen bee is represented as the transitions in a state 
space (search space), with the queen probabilistically mat¬ 
ing with the drone encountered at each state. The prob¬ 
ability of mating is determined by the speed and energy 
of the queen, and the fitness of the drone. Karaboga and 
Basturk proposed an Artificial Bee Golony optimization 
(ABG) in [TB] . ABG classifies the bees in a hive into three 
types: “scout bees” that randomly fly without guidance, 
“employed bees” that search the neighborhood of their po¬ 
sitions, and “onlooker bees” that use the population fitness 
to select a guiding solution for exploitation. The algo¬ 
rithm balances exploration and exploitation by means of 
using employed and onlooker bees for local search, and the 
scout bees for global search. It also demonstrates satisfac¬ 
tory performance in applications [19] [20] . 

Besides the bees, other organisms have also been widely 
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studied [5]. Krishnanand and Ghose proposed a Glow¬ 
worm Swarm Optimization (GSO) [H] based on the be¬ 
havior of the firefly. In GSO, each firefly randomly se¬ 
lects a neighbor according to its luminescence and moves 
toward it. In general the fireflies are more likely to get 
interested in others that glow brighter. As the movement 
is only conducted locally using selective neighbor informa¬ 
tion, the firefly swarm is able to divide into disjoint sub¬ 
groups to explore multiple optima. Another firefly-based 
technique is proposed by Yang et al. |22j . He reformu¬ 
lated the co-movement pattern of fireflies and employed it 
in optimization. Passino devised a Bacterial Foraging Op¬ 
timization (BFO) [53] based on the bacterial chemotaxis. 
In BFO, possible solutions to the optimization problem are 
represented by a colony of bacteria. It consists of three 
schemes, i.e., chemotaxis, reproduction, and elimination- 
dispersal. The exploitation task is performed using the 
first two schemes and the last one contributes to explo¬ 
ration [53]. Researchers have also devised swarm intelli¬ 
gence algorithms based on other organisms and they can 
also generate satisfactory optimization performance |3]. 
To the best of our knowledge, only one spider-inspired 
metaheuristic aiming at solving optimization problem has 
been proposed, i.e. the Social Spider Optimization |5S| 
devised by Guevas et al., which divides the spiders into 
different genders and mimics the mating behavior for op¬ 
timization. However, our proposed algorithm is totally dif¬ 
ferent from this algorithm in their biological backgrounds, 
motivations, implementations, and search behaviors. We 
will further reveal the differences in Section [SHj 

3. Social Spider Algorithm 

In SSA, we formulate the search space of the optimiza¬ 
tion problem as a hyper-dimensional spider web. Each po¬ 
sition on the web represents a feasible solution to the opti¬ 
mization problem and all feasible solutions to the problem 
have corresponding positions on this web. The web also 
serves as the transmission media of the vibrations gener¬ 
ated by the spiders. Each spider on the web holds a posi¬ 
tion and the quality (or fitness) of the solution is based on 
the objective function, and represented by the potential 
of finding a food source at the position. The spiders can 
move freely on the web. However, they can not leave the 
web as the positions off the web represent infeasible solu¬ 
tions to the optimization problem. When a spider moves 
to a new position, it generates a vibration which is propa¬ 
gated over the web. Each vibration holds the information 
of one spider and other spiders can get the information 
upon receiving the vibration. 

3.1. Spider 

The spiders are the agents of SSA to perform opti¬ 
mization. At the beginning of the algorithm, a pre-defined 
number of spiders are put on the web. Each spider s holds 
a memory, storing the following individual information: 


• The position of s on the web. 

• The fitness of the current position of s. 

• The target vibration of s in the previous iteration. 

• The number of iterations since s has last changed its 
target vibration. 

• The movement that s performed in the previous it¬ 
eration. 

• The dimension masl|^that s employed to guide move¬ 
ment in the previous iteration. 

The first two types of information describe the individual 
situation of s, while all others are involved in directing s 
to new positions. The detailed scheme of movement will 
be elaborated in Section 1331 

Based on observations, spiders are found to have very 
accurate senses of vibration. Furthermore, they can sepa¬ 
rate different vibrations propagated on the same web and 
sense their respective intensities m- In SSA, a spider will 
generate a vibration when it reaches a new position differ¬ 
ent from the previous one. The intensity of the vibration is 
correlated with the fitness of the position. The vibration 
will propagate over the web and other spiders can sense 
it. In such a way, the spiders on the same web share their 
personal information with others to form a collective social 
knowledge. 

3.2. Vibration 

Vibration is a very important concept in SSA. It is one 
of the main characteristics that distinguish SSA from other 
metaheuristics. In SSA, we use two properties to define 
a vibration, namely, the source position and the source 
intensity of the vibration. The source position is defined 
by the search space of the optimization problem, and we 
define the intensity of a vibration in the range [0,-l-oo). 
Whenever a spider moves to a new position, it generates 
a vibration at its current position. We define the position 
of spider a at time t as Pa{t), or simply as Pa if the time 
argument is t. We further use I{Pa, Pb,t) to represent 
the vibration intensity sensed by a spider at position P^ 
at time t and the source of the vibration is at position 
Pa. With these notations we can thus use I{Ps, Ps,t) to 
represent the intensity of the vibration generated by spider 
s at the source position. This vibration intensity at the 
source position is correlated with the fitness of its position 
/(Pg), and we define the intensity value as follows: 

I{Ps,Ps,t) = \og{j^py^ + l) (1) 

where C is a confidently small constant such that all pos¬ 
sible fitness values are larger than C. Please note that we 
consider minimization problems in this paper. The design 
of Q takes the following issues into consideration: 


^The dimension mask is a 0-1 binary vector of length D, where D 
is the dimension of the optimization problem 
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• All possible vibration intensities of the optimization 
problem are positive. 

• The positions with better fitness values, i.e. smaller 
values for minimization problems, have larger vibra¬ 
tion intensities than those with worse fitness values. 

• When a solution approaches the global optimum, the 
vibration intensity would not increase excessively, 
and cause malfunctioning of the vibration attenu¬ 
ation scheme. 

As a form of energy, vibration attenuates over distance. 
This physical phenomenon is accounted for in the design 
of SSA. We define the distance between spider a and b as 
D{Pa, Pb) and we use 1-norm (Manhattan distance) to 
calculate the distance, i.e., 

D{Pa,Pb)=\\Pa-Pb\\l- (2) 

The standard deviation of all spider positions along each 
dimension is represented by cr. With these definitions, we 
further define the vibration attenuation over distance as 
follows: 


IiPa,Pb,t) = I{Pa,Pa,t) X exp(-. (3) 

a XTa 

In the above formula we introduce a user-controlled pa¬ 
rameter Tfl S (0,oo). This parameter controls the atten¬ 
uation rate of the vibration intensity over distance. The 
larger is, the weaker the attenuation imposed on the 
vibration. 

3.3. Search Pattern 

Here we demonstrate the above ideas in terms of an 
algorithm. There are three phases in SSA: initialization, 
iteration, and final. These three phases are executed se¬ 
quentially. In each run of SSA, we start with the initializa¬ 
tion phase, then perform searching in an iterative manner, 
and finally terminate the algorithm and output the solu¬ 
tions found. 

In the initialization phase, the algorithm defines the 
objective function and its solution space. The value for 
the parameter used in SSA is also assigned. After setting 
the values, the algorithm proceeds to create an initial pop¬ 
ulation of spiders for optimization. As the total number of 
spiders remains unchanged during the simulation of SSA, 
a fixed size memory is allocated to store their informa¬ 
tion. The positions of spiders are randomly generated in 
the search space, with their fitness values calculated and 
stored. The initial target vibration of each spider in the 
population is set at its current position, and the vibration 
intensity is zero. All other attributes stored by each spider 
are also initialized with zeros. This finishes the initializa¬ 
tion phase and the algorithm starts the iteration phase, 
which performs the search with the artificial spiders cre¬ 
ated. 


In the iteration phase, a number of iterations are per¬ 
formed by the algorithm. In each iteration, all spiders on 
the web move to a new position and evaluate their fitness 
values. Each iteration can be further divided into the fol¬ 
lowing sub-steps: fitness evaluation, vibration generation, 
mask changing, random walk, and constraint handling. 

The algorithm first calculates the fitness values of all 
the artificial spiders on different positions on the web, and 
update the global optimum value if possible. The fitness 
values are evaluated once for each spider during each it¬ 
eration. Then these spiders generate vibrations at their 
positions using 0 . After all the vibrations are generated, 
the algorithm simulates the propagation process of these 
vibrations using ([^. In this process, each spider s will re¬ 
ceive \pop\ different vibrations generated by other spiders 
where pop is the spider population. The received infor¬ 
mation of these vibrations include the source position of 
the vibration and its attenuated intensity. We use V to 
represent these \pop\ vibrations. Upon the receipt of V, s 
will select the strongest vibration from V and com¬ 
pare its intensity with the intensity of the target vibration 
stored in its memory, s will store as if the 
intensity of is larger, and c^, or the number of it¬ 
erations since s has last changed its target vibration, is 
reset to zero; otherwise, the original Vtar is retained and 
Cs is incremented by one. We use P* and to repre¬ 
sent the source positions of V and Vtar, respectively, and 
i = {1,2, • • • , \pop\}. 

The algorithm then manipulates s to perform a random 
walk towards 1 ;*“’". Here we utilize a dimension mask to 
guide the movement. Each spider holds a dimension mask 
m, which is a O-I binary vector of length D and D is the 
dimension of the optimization problem. Initially all values 
in the mask are zero. In each iteration, spider s has a 
probability of 1 —Pc®° to change its mask where Pc S (0,1) 
is a user-defined attribute that describes the probability of 
changing mask. If the mask is decided to be changed, each 
bit of the vector has a probability of Pm to be assigned with 
a one, and I — Pm to be a zero. Pm is also a user-controlled 
parameter defined in (0,1). Each bit of a mask is changed 
independently and does not have any correlation with the 
previous mask. In case all bits are zeros, one random value 
of the mask is changed to one. Similarly, one random bit 
is assigned to zero if all values are ones. 

After the dimension mask is determined, a new follow¬ 
ing position P{° is generated based on the mask for s. 
The value of i-th dimension of the following position P{°j 
is generated as follows. 




ms,I = 0 

ms,i = I ’ 


(4) 


where r is a random integer value generated in [1, |pop|], 
and ms,i stands for the i-th dimension of the dimension 
mask m of spider s. Here the random number r for two 
different dimensions with ms,i = 1 is generated indepen¬ 
dently. 


4 



With the generated P{°, s performs a random walk to 
the position. This random walk is conducted using the 
following equation. 

Psit+1) = Ps + iPs-Ps{t-l))xr+{P{°-Ps)QR, (5) 


where © denotes element-wise multiplication and i? is a 
vector of random float-point numbers generated from zero 
to one uniformly. Before following P{°, s first moves along 
its previous direction, which is the direction of movement 
in the previous iteration. The distance along this direc¬ 
tion is a random portion of the previous movement. Then 
s approaches P{° along each dimension with random fac¬ 
tors generated in (0,1). This random factor for different 
dimensions are generated independently. After this ran¬ 
dom walk, s stores its movement in the current iteration 
for the next iteration. This ends the random walk sub¬ 
step. 

The final sub-step of the iteration phase is the con¬ 
straint handling. The spiders may move out of the web 
during the random walk step, which causes the constraints 
of the optimization problem to be violated. There are 
many methods to handle the boundary constraints in the 
previous literature, and the random approach, absorbing 
approach, and the reflecting approach are three most widely- 
adopted methods [5S]. In this paper we adopt the re¬ 
flecting approach for constraint handling and produce a 
boundary-constraint-free position Ps(t -|- 1) by 


Ps,i{t + 1) 


(Xi Ps^i) XT if Psji(t T 1) © Xi , , 
(Ps,j - X r if Ps,i(t -I- 1) < ^ ’ 


where xi is the upper bound of the search space in the i-th 
dimension, and Xi is the lower bound of the corresponding 
dimension, r is a random floating point number generated 
in (0,1) 

The iteration phase loops until the stopping criteria is 
matched. The stopping criteria can be defined as the max¬ 
imum iteration number reached, the maximum CPU time 
used, the error rate reached, the maximum number of it¬ 
erations with no improvement on the best htness value, or 
any other appropriate criteria. After the iteration phase, 
the algorithm outputs the best solution with the best ht¬ 
ness found. The above three phases constitute the com¬ 
plete algorithm of SSA and its pseudo-code can be found 
in Algorithm 


Algorithm 1 Social Spider Algorithm 
1: Assign values to the parameters of SSA. 

2: Create the population of spiders pop and assign mem¬ 
ory for them. 

3: Initialize for each spider. 

4: while stopping criteria not met do 
5: for each spider s in pop do 

6: Evaluate the htness value of s. 

7: Generate a vibration at the position of s. 

8 : end for 

9: for each spider s in pop do 

10: Calculate the intensity of the vibrations V 

generated by all spiders. 

11: Select the strongest vibration from V. 

12: if The intensity of is larger than then 

13: Store as 

14: end if 

15: Update Cs- 

16: Generate a random number r from [0,1). 

17: if r > then 

18: Update the dimension mask m-g. 

19: end if 

20: Generate P{°. 

21: Perform a random walk. 

22: Address any violated constraints. 

23: end for 

24: end while 

25: Output the best solution found. 


3.4- Differences between SSA and Other Evolutionary Com¬ 
putation Algorithms 

A number of swarm intelligence algorithms have been 
proposed in the past few decades. Among them, PSO 
and AGO are the two most widely employed and studied 
algorithms. SSA may also be classihed as a swarm intel¬ 
ligence algorithm, but it has many differences from PSO 
and AGO, elaborated below. 

PSO, like SSA, was originally proposed for solving con¬ 
tinuous optimization problems. It was also inspired by 
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animal behavior. However, the first crucial difference be¬ 
tween SSA and PSO is in individual following patterns. 
In PSO, all particles follow a common global best position 
and their own personal best position. However in SSA, 
all spiders follow positions constructed by others’ current 
positions and their own historical positions. These fol¬ 
lowing positions are not guaranteed to be visited by the 
population before, and different spiders can have different 
following positions. Since the global best position and spi¬ 
ders’ current positions differ greatly during most time of 
the optimization process, these two following patterns lead 
to different searching behaviors. This may weaken the con¬ 
vergence ability of SSA but can potentially strengthen the 
capability of solving multi-modal optimization problems 
with a great number of local optimums. 

Besides the difference in the following pattern, the dif¬ 
ference in their biology backgrounds is also very signifi¬ 
cant. PSO was designed based on the model of coordi¬ 
nated group animal motions of flocks of birds or schools 
of fishes. This model serves as the design metaphor of 
PSO. SSA is inspired by the social spider foraging strat¬ 
egy, which belongs to the scope of general social animal 
searching behavior. We use a general IS model as the 
design framework. This difference is also a major distin¬ 
guishing feature of SSA from other proposed algorithms. 

A third difference between SSA and the original formu¬ 
lation of PSO is in the information propagation method. 
In PSO, the information propagation method is neglected, 
and each particle is assumed to be aware of all the infor¬ 
mation of the system without loss. Although the informa¬ 
tion validity ranges are considered in some recent variants 
of PSO, the information loss characteristic is still a unique 
feature that distinguishes SSA from PSO variants. In SSA 
we model the information propagation process through the 
vibrations on the spider web. This process forms a gen¬ 
eral knowledge system with information loss. Although 
there is still no research on how the information loss will 
impact the social foraging strategy employed in optimiza¬ 
tion, it is possible that this information loss system par¬ 
tially contributes to the performance improvement of SSA 
over PSO. 

Another difference is that in PSO, the common knowl¬ 
edge of the group is all about the best particle in the sys¬ 
tem. All remaining particles in the system do not consti¬ 
tute the shared information of the group, which may lead 
to neglecting some valuable information of the population. 
In SSA, each spider generates a new piece of information 
and propagates the information to the whole population. 
There are also some PSO variants that shares personal best 
position information with the population, but the main fo¬ 
cus of PSO is on the best positions of the individuals and 
the population. In SSA, the information generated and 
propagated with the vibrations are the current positions 
instead of the best-in-history positions, which may differ 
greatly with the search. 

Although both SSA and AGO draw their inspirations 
from the social animal foraging strategy, there are still 


some obvious differences. The foraging frameworks adopted 
by the two algorithms are quite different: AGO utilizes the 
ant foraging behavior to perform optimization. Ants find 
food by laying down pheromone trails and collectively es¬ 
tablishing positive feedbacks which bias the later path se¬ 
lection, while spiders sense the vibration propagated by 
the spider web to locate the prey. Another difference is 
the presentation of feasible solutions. In SSA we use the 
positions on the spider web to represent feasible solutions. 
Similar representations have also been widely adopted in 
the swarm intelligence algorithm. Meanwhile, AGO uses 
the path between the ant hive and food sources to rep¬ 
resent solutions to the optimization problems. Addition¬ 
ally, AGO was originally designed to solve combinatorial 
problems. Although in recent years there are AGO-variant 
algorithms designed mainly to solve continuous problems 
P7] , the performance is not as good as the performance of 
the original AGO in solving combinatorial problems like 
the Traveling Salesman Problem. There are also infor¬ 
mation propagation and searching pattern differences be¬ 
tween SSA and AGO as described above. 

There are also some other swarm intelligence algorithms 
proposed to solve continuous problems, and SSA has some 
unique characteristics. In most swarm intelligence algo¬ 
rithms, e.g., ABG and GSO, the populations are struc¬ 
tured into different types. Different types of individuals 
perform different jobs and the whole population cooperates 
to search the solution space. However in SSA, all individ¬ 
uals (spiders) are equal. Each performs all the tasks that 
would be executed by multiple types of the populations 
in other algorithms. If we put SSA into the conventional 
framework, it has the feature that the different types of 
individuals can transform into other types very smoothly 
and without the guidance of the user, which may poten¬ 
tially contribute to the performance improvement. 

As to Social Spider Optimization (SSO) [2S], differ¬ 
ences lie in all aspects of the algorithm design. A most 
important difference is that in SSO the spiders are clas¬ 
sified by gender. Male and female spiders have different 
searching operations. However, the spiders in SSA share 
the same searching operation, significantly reducing the 
effort in implementation. SSA also incorporates the infor¬ 
mation propagation model into its algorithm design, and 
thus the social spider population in SSA fits the IS model. 
Besides, SSA imitates the foraging behavior of social spi¬ 
ders, while SSO imitates the mating behavior of social 
spiders. The differences in algorithm implementation are 
more patent. In SSO there are three spider movement 
operators executed first in parallel and then in sequence. 
The moving pattern of the third operator highly depends 
on the first two operators. This design may potentially 
increase the difficulty of analyzing the search behavior of 
the algorithm. In SSA we implement one random move op¬ 
erator, which combines both exploration and exploitation 
behaviors in one move. In our design, the search behav¬ 
ior is controlled by the parameters, thus providing a clear 
view on the search behavior of the algorithm. The impact 
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of different parameters on the optimization performance 
of SSA is further illustrated in Section [331 

Although EAs, like GA and ES, are also population- 
based algorithms, and inevitably share some similarities 
with the population-based SSA, they are quite different 
general-purpose metaheuristics. They are inspired by com¬ 
pletely different biological disciplines. EAs usually employ 
different recombination and decomposition operators to 
manipulate the solutions, which imitate the regeneration 
of an organism. 

As stated above, although we still do not know the ex¬ 
act impact of information loss on the optimization process, 
this feature of SSA may contribute to the optimum search 
in some complex multimodal optimization problems. The 
uniform structure of the population is another potential 
advantage of SSA. In addition, the unique searching pat¬ 
tern and its underlying social animal foraging strategy as 
well as the IS foraging model contribute to the overall per¬ 
formance of SSA. 

3.5. Adjusting SSA Parameters 

Choosing proper parameters of SSA for numerical and 
real-world optimization problems can be time-consuming. 
The trial-and-error scheme, or a parameter sweep test, 
may reveal the best achievable performance over the pa¬ 
rameter space at the expense of high computational cost. 
In real-world optimization problems, evaluating the fitness 
function may take a long time, much longer than evalu¬ 
ating our benchmark functions, and one evaluation may 
take several seconds or even minutes to finish, rendering 
trial-and-error schemes impractical for parameter tuning. 
As alternatives, researchers have proposed some schemes 
to replace the trial-and-error parameter selection scheme. 
These schemes can generally be classified into three groups 

m- 

• Fixed parameter schemes select a parameter combi¬ 
nation before the simulation using empirical or theo¬ 
retical knowledge of the characteristics of the param¬ 
eters. This combination remains constant through¬ 
out the whole search [29] [30| . 

• Deterministic parameter schemes use some pre-defined 
rules to change the parameter values throughout the 
search pU]|5T]. 

• Adaptive parameter schemes change the parameter 
values by adaptively learning the impact of changing 
parameters on the searching performance through¬ 
out the search |^. Some schemes encode the pa¬ 
rameters into the solution and evolve the parameters 
together with the population |33|. 

In this paper we use the fixed parameter scheme to test 
the performance of SSA compared with other algorithms. 
We also use this scheme to perform a preliminary param¬ 
eter sensitivity analysis in order to deduce some rules of 
thumb on choosing parameters that can consistently lead 


to satisfactory results on a wide range of functions with 
different characteristics. This test can also discover some 
of the features of the parameters when solving different 
kinds of optimization problems. We carry out extensive 
simulations on our benchmark functions, which cover a 
wide range of optimization problems. Thus, the derived 
rules of thumb can be expected to give generally good 
performance on unknown problems. 

In SSA we employ three user-controlled parameters to 
guide the searching behaviour, namely, 

• ra'. This parameter defines the rate of vibration at¬ 
tenuation when propagating over the spider web. 

• Pc'. This parameter controls the probability of the 
spiders changing their dimension mask in the random 
walk step. 

• Pm'- This parameter defines the probability of each 
value in a dimension mask to be one. 


In this section, we employ five 10-dimensional bench¬ 
mark functions to investigate the impact of these param¬ 
eters on the performance of SSA. These functions are the 
Sphere, Schwefel 2.22, Rastrigin, Ackley, and Griewank 
functions and the detailed definition of these functions can 
be found in Section |4| All benchmark functions are not 
shifted. The value of is selected from the set {|^|^, 1,2,3, 
and the values of Pc and Pm are both selected from the 
set {0.01,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,0.99}. So for 
each function/parameter pair there are 11 data points for 
analysis. The stopping criteria is set to 100 000 evaluations 
and each function is tested for 20 times. The population 
size is set to 10 in accordance with our later simulation/ 

The mean results are plotted in Eig. [^with dots, and the 
second-order polynomial regression curve for each func¬ 
tion is also plotted for demonstration. As the Sphere and 
Schwefel 2.22 functions are uni-modal functions while the 
other three are multi-modal functions, we can obtain some 
interesting observations from the mean results and the re¬ 
gression curve. 

Erom Fig. [^we can see SSA is very robust in solving 
unimodal problems when the three parameters change. All 
six sub-figures indicate outstanding performance in terms 
of fitness value. A general conclusion is that the value of 
Pc shall not be set to a very large value, e.g. 0.99. 

When considering the multimodal problems, we can 
observe some obvious impact of the parameters on the 
performance. In terms of the vibration attenuation rate 
Tq, both the mean results and the regression curves favor 
relatively small values. A generally preferred value is one, 
while random values selected from (0, 3] shall also be able 
to yield good results. 

While ra presents relatively stable results in the mul¬ 
timodal problem tests, Pc has different impacts on differ¬ 
ent problems. Rastrigin function favors a small Pc while 
Rosenbrock function prefers a medium-large one. It seems 
that most Pc values can generate good results on Ackley 
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Figure 1: Parameter test results on Pc, and pm- 


function. So the selection of Pc highly depends on the 
nature of the optimization problems to be solved. 

As to Pm, the simulation results reveal a similar ten¬ 
dency with Ta that the multimodal problems seem to prefer 
a relatively small Pm somewhere near 0.1. So we adopt the 
parameter combination = 1, pc = 0.7, and pm — 0.1 in 
all of our later simulations. Please note that this param¬ 
eter combination is not guaranteed to be the best one for 
solving all optimization problems, and parameter tuning 
is essential to address unfamiliar problems. 

Please note that the parameter sensitivity analysis in 
this paper is a preliminary one. In this test only one pa¬ 
rameter is tested while the remaining two is set unchanged, 
i.e., Ta = 1, Pc = 0.7, and Pm = 0.1. Although it is 
not guaranteed that this set of parameters works best for 
the benchmark functions, it is one that yields outstanding 
performance, which will be demonstrated in Section]^ A 
complete parameter sensitivity analysis is one of the future 
research topics of SSA. 

4. Benchmark Problems and Evaluation Method 

In order to benchmark the performance of SSA, we 
conduct simulations on 25 different benchmark functions. 
These benchmark functions are all the base functions from 
the latest Competition on Real-Parameter Single Objec¬ 
tive Optimization Problems at CEC 2013 [34] and CEC 


2014 |3S]. The benchmark functions can be classified into 
four groups: 

• Group I: / 1-/5 are unimodal functions. 

• Group II: fe-fi 5 are multimodal functions. 

• Group III: /i 6“/20 are rotated multimodal functions 
whose base functions belong to Group II functions. 

• Group IV: / 21-/25 are hybrid multimodal functions 
whose base functions belong to Group I - III func¬ 
tions. 

The benchmark functions are listed in Table [TJ All bench¬ 
mark functions, except /13 Schwefel’s Problem 2.26, are 
shifted minimization problems and the search ranges are 
scaled to [— 100 , 100 ]”, where n is the dimension of the 
problem. Group I functions are used to test the fast¬ 
converging performance of SSA. Group II functions all 
have a large number of local minima points, and can be 
used to test the ability of SSA to jump out of local optima 
and avoid pre-mature convergence. Group III functions 
are more complex than other functions and can push the 
searching capability of SSA to a limit. Group IV func¬ 
tions are employed to test the optimization performance of 
handling problems consisting of different subcomponents 
with different properties. The detailed implementation of 
Group IV functions, i.e. hybrid multimodal functions, can 
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Table 1: Benchmark Functions 
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* o is a shifting vector and AT is a transformation matrix, o and AT can be obtained from m- 
■*" rand{) is a random number uniformly generated in (0,1). 


be found in [3S]. In Table[^we only list the subcomponents 
of the hybrid functions. 

All benchmark functions locate their global minimum 
values at zero, and the fitness values smaller than 10“® are 
considered as 10“® as required in [33] and |3S|. We test 
the benchmark functions in 10, 30, and 50 dimensions to 
draw empirical conclusion on the scalability of SSA, and 
each function is tested for 51 runs [35|. In each run, we 
use a maximum number of 10^ x n function evaluations as 
the termination criteria. In order to meet the requirement 
set by |35j , we use one fixed combination of parameters for 
SSA in the simulation of all groups of functions. The pop¬ 
ulation size is n, and other parameters are set according 
to the analysis in Section [33] 

To evaluate the performance of SSA, we compare the 
simulation results with the state-of-the-art algorithms in 
solving real-parameter optimization problems, including 
the variances of Co-variance Matrix Adaptation Evolution 
Strategies (CMA-ES) |36|, adaptive Differential Evolution 
algorithms (DE) |37|, and a Global and Local real-coded 
Genetic Algorithm (GL-25) |3H|. CMA-ES and DE vari¬ 
ants are arguably the most successful optimization algo¬ 


rithms current in use |3Hj . In the latest Competition on 
Real-Parameter Single Objective Optimization Problems 
at CEC 2013, their variant algorithms possess nine po¬ 
sitions in the top ten best performing algorithms. For 
DE variants, we select JADE |3S| and SaDE [33] for com¬ 
parison due to their excellent performance demonstrated 
in the CEC 2013 competition. All source codes are ob¬ 
tained from the original author. We make some minor 
changes to adapt them to our benchmark functions, but 
the main body and logic of the algorithms are untouched. 
The stopping criteria for all the compared algorithms are 
set to 10"^ X n function evaluations, and the parameters of 
these algorithms are set according to the recommendation 
made in the corresponding literature, i.e., CMA-ES in [IT] . 
JADE in [301, SaDE in |53|, and GL-25 in [35] . 

Besides these state-of-the-art algorithms, we also per¬ 
formed simulations with other famous algorithms, namely, 
Real-Coded Genetic Algorithm [32], Adaptive Particle Swarm 
Optimization [33] , Artificial Bee Colony Optimization [18] , 
Firefly Algorithm [33], Cuckoo Search [33], and Group 
Search Optimizer |45| . As their overall performance in 
terms of best fitness values achieved and convergence speed 
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is not comparable to SSA and the other four algorithms, 
the detailed simulation results will not be presented in this 
paper. 

5. Numerical Experiments and Results 

In this section we present the simulation results of SSA 
on the benchmark functions identified in Section m We 
perform comparison among SSA and other algorithms and 
give statistical analysis on the simulation results. 

5.1. Experimental Comparison with Other State-of-the-Art 
Algorithms 

We first conduct a series of simulations on the 30- 
dimension optimization problems using SSA and other state- 
of-the-art algorithms. The simulation results are plotted 
in Fig. 1^ and Table reports the means and standard 
deviations of the optimal fitness values achieved, with the 
best mean result shaded. Besides, we also perform a series 
of Wilcoxon rank sum tests on the null hypothesis that 
SSA performs similarly with other algorithms when solv¬ 
ing each benchmark function. The statistical test result 
at 95% significance level is presented next to the standard 
deviation of the corresponding algorithm, where a © in¬ 
dicates that SSA performs significantly better than the 
tested algorithm on the specified function, a © indicates 
that SSA performs not as good as the tested algorithm, 
and a © means that the Wilcoxon rank sum test cannot 
distinguish between the simulation results of SSA and the 
tested algorithm. The counts of the benchmark functions 
that fall in these situations are shown at the bottom of the 
table. 

From Table we observe that: 

• SSA generally outperforms all compared algorithms 
in terms of the statistical test. Among all 25 func¬ 
tions, SSA generates better simulation results in 16, 
12, 14, and 13 functions compared with CMA-ES, 
JADE, SaDE, and GL-25, respectively. If we take 
those functions with similar results, the advantage is 
more obvious: SSA performs no worse than CMA- 
ES, JADE, SaDE, and GL-25 in 22, 19, 24, and 22 
functions, respectively. 

• In the first group of benchmark functions, all com¬ 
pared algorithms can obtain the global optimum val¬ 
ues of / 1-/4 in all runs, which means that the final 
result test cannot reveal the best-performing algo¬ 
rithms. We shall further employ the convergence test 
to analyse the performance of compared algorithms. 

• The performance of SSA in solving Group II multi¬ 
modal optimization problems is superior, and it gen¬ 
erates 7 best mean results out of the total 10. The 
numbers for CMA-ES, JADE, SaDE, and GL-25 are 
1, 4, 1, and 4, respectively. Besides, the mean re¬ 
sults of the three functions where SSA is not the 


best performing one are still very competitive and 
comparable to the best results. 

• SSA is not as competitive in solving rotated mul¬ 
timodal functions as it does in unimodal and mul¬ 
timodal functions. However, a careful investigation 
on the mean results shows that the performance of 
SSA is still comparable to all compared algorithm. 
The reason of this phenomenon may be that dur¬ 
ing the searching process of SSA, no correlation ma¬ 
trix or differential vectors are employed to assist the 
blind search in the solution space as in CMA-ES and 
JADE. However, this disadvantage can be overcome 
by employing these mentioned schemes into SSA or 
via hybrid algorithms. This is a potential future re¬ 
search direction. 

• SSA is very powerful at solving hybrid functions where 
different dimensions of the objective functions can be 
un-related and there is no additional information like 
correlation matrix available. In Group IV tests, SSA 
achieved four out of the total five best mean results. 

5.2. Scalability Test 

In addition to the 30-dimension benchmark function 
tests, we also performed a series of simulations on both 
10- and 50-dimension benchmarks to test the scalability of 
SSA. To make a thorough comparison, we also employed 
the compared algorithms in this test. The simulation re¬ 
sults are presented in Tables and using the same for¬ 
mat and symbols as in Table From the results we have 
the following observations: 

• The advantage of SSA over compared algorithms are 
confirmed. SSA achieved all the best mean results in 
10-D Group H tests and 50-D Group IV tests, and 
have satisfactory performance compared with other 
algorithms in all other groups. 

• This advantage is also supported by the statistical 
test. SSA can generate better results than CMA-ES, 
JADE, SaDE, and GL-25 in 18, 15, 19, and 8 10-D 
functions, respectively. The corresponding numbers 
for 50-D functions are 17, 9, 12, and 15, respectively. 
The statistical results will favor SSA more if we also 
take those tests that have similar performance into 
account. 

• From the simulation results we can see GL-25 per¬ 
forms very well in 10-D hybrid functions, JADE per¬ 
forms very well in 50-D multimodal functions, when 
compared with the remaining three algorithms. How¬ 
ever, SSA can always outperform them in these tests, 
which indicate the superior scalability of SSA. 

5.3. Convergence Test 

As stated in Section |5.11 the final result comparison 
cannot completely describe the searching performance of 
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Figure 2: Box plot of raw simulation results. 
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Figure 3: Median convergence test results. 
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an algorithm. So we further conduct a convergence test 
on the five compared algorithms on each 30-D benchmark 
function. We employ the raw simulation data generated in 
Section [Sd] As each function is tested for 51 runs for each 
algorithm, we select the convergence data of the run which 
generates the median final result. The convergence data of 
the five compared algorithms are plotted in Fig. The x- 
axis is the function evaluations consumed, and the y-axis 
is the best-so-far fitness values found. The convergence 
plots lead to the following observations: 

• The convergence speed of SSA in solving unimodal 
optimization problems is not as fast as CMA-ES, 
JADE, and SaDE. This is because SSA performs 
exploitation and exploration simultaneously during 
the random walk process, and the former contributes 
to a fast convergence speed but the latter obstructs 
the population from moving into a small region in 
the search space. However, the degree of exploration 
in the random walk process can be effectively con¬ 
trolled by the algorithm parameters as revealed in 
Section |3.5[ So it is highly possible that SSA can 
also achieve a comparable convergence speed with 
other algorithms if we adopt a set of suitable pa¬ 
rameters designed for solving unimodal optimization 
problems. 

• When solving multi-modal optimization problems, 
SSA generally converges as fast as or even faster than 
the compared algorithms. This phenomenon can be 
clearly observed in the Group II tests, where SSA 
mostly generates a similar convergence curve with 
others. 

• The advantage of combining exploitation and explo¬ 
ration in one searching process is revealed in the con¬ 
vergence plot. Take /e as an example, almost all al¬ 
gorithms (except JADE) are trapped in local optima 
shortly after the start of searching. However, SSA 
managed to jump out of the optimum and success¬ 
fully found the global optimum by the end of search¬ 
ing. The manifestation of this searching characteris¬ 
tic is that the algorithm convergences relatively slow 
at the beginning, but then very fast as the search 
continues. This phenomenon can also be observed in 
some other instances, e.g. / 13 , / 22 , and / 24 . 

5.4- Reliability Test 

Another test which can examine the reliability of stochas¬ 
tic algorithms is the success rate test and this test has 
been adopted in many previous work, e.g. [32] [13]. In 
this part, we aim at further comparing and visualizing 
the performance of all compared algorithms on all tested 
benchmark functions by plotting the empirical cumulative 
distribution of success rates. The plots are presented in 
Fig. i where the x-axis is different success thresholds, 
and the y-axis is the overall success rate. A simulation 
run is considered successful if and only if the best-found 


fitness value is smaller than or equal to the “successful 
threshold”. The overall success rate is calculated by the 
number of successful runs under a specific success thresh¬ 
old divided by the total number of runs. Thus a larger 
overall success rate implies a more reliable algorithm. By 
comparing the distribution curves of different algorithms, 
we can have a general view on which algorithm is reliable 
at solving general optimization problems. 

From the plots we have the following observations: 

• SSA is generally more reliable than the compared al¬ 
gorithms. The advantage is very significant in 10-D 
and 50-D benchmark functions. While the reliabil¬ 
ity of SSA is very similar to JADE in small suc¬ 
cess threshold regions, SSA regains the lead in large 
threshold regions. 

• The plot also demonstrates the convergence charac¬ 
teristic of SSA. SSA is able to find the global op¬ 
timum whenever the algorithm is able to locate a 
relatively small region near it. This conclusion is 
drawn based on the observation that almost all runs 
that obtained a fitness value smaller than 10 “"* in 
30-D functions, and 10“^ in 50-D functions are able 
to converge to the global optimum point at 10 “®. 

5.5. Computational Complexity 

Besides the previous tests, the computational complex¬ 
ity is also a major factor for evaluating the efficiency of 
an evolutionary computation algorithm. In this paper, 
we employ the method stated in [33] and [35] to analyze 
the computational complexity of the compared algorithms. 
We use fe in Table [l] which is the major component of the 
testing methodology suggested by [35], as the benchmark 
evaluation function, and the complexity analysis result is 
as follows. The complexity values of CMA-ES, JADE, 
SaDE, GL-25, and SSA are 40.47, 34.69, 75.72, 63.51, and 
44.18, respectively. A smaller complexity value means that 
the algorithm is less complex, which leads to a relatively 
faster speed in execution under the same condition. From 
the results we can see that although SSA is slightly more 
complicated than JADE and CMA-ES, their complexities 
are comparable. In addition, these three algorithms are 
significantly less computationally complex than SaDE and 
GL-25. 

5.6. Discussion 

As stated by the No-Free-Lunch (NFL) Theorem [33] , 
all meta-heuristics that search for extrema shall perform 
exactly the same when all possible objective functions are 
evaluated and averaged. It is further elaborated that it 
is theoretically impossible to have a best general-purpose 
universal optimization technique m- Superior performing 
algorithms are available if particular classes of problems 
are considered m, or general but real-world ones |48j . 


13 


50-D 





Figure 4: Empirical cumulative distribution of success rate test results. 


However, the total nuniber of possible problems are so 
huge that there is still much room to develop new algo¬ 
rithms. Though existing meta-heuristics have great suc¬ 
cess in solving many optimization problems, it is always 
worthwhile to develop new searching methodologies with 
superior performance in particular classes of problems |49j . 
This is the motivation for us to propose SSA for solving 
global numerical optimization problems. 

6. Conclusion 

In this paper we proposed a novel social spider algo¬ 
rithm to solve global optimization problems. This algo¬ 
rithm is based on the foraging behavior of social spiders 
and the information-sharing foraging strategy. SSA is con¬ 
ceptually simple and relatively easy to implement. SSA 
can tackle a wide range of different continuous optimiza¬ 
tion problems and has the potential to be employed to 
solve real-world problems. 

In order to evaluate the performance of SSA, we adopted 
a set of 25 benchmark functions which cover a large vari¬ 
ety of different optimization problem types. We compared 
SSA with the state-of-the-art optimization algorithms, namely, 
CMA-ES, JADE, SaDE, and GL-25. These algorithms 
have been employed to solve a large set of different bench¬ 
mark optimization functions and real-world problems, and 
demonstrated outstanding performance. The results show 
that the performance of SSA is outstanding compared with 
the above listed algorithms in all three different groups of 
functions. This conclusion was supported by both the sim¬ 
ulation results and the statistics of the simulation data. 

Euture research on SSA can be divided into three cat¬ 
egories: scheme research, algorithm research, and real- 
world application. The random walk scheme in the current 
SSA may be further improved using advanced optimiza¬ 
tion techniques and hybrid algorithms with deterministic 
heuristics or local search algorithms. New schemes can 
also be applied in the searching process of SSA for per¬ 
formance improvement. In terms of algorithm research. 


SSA has the potential to be applied to solve combinato¬ 
rial problems. We note that some other swarm intelli¬ 
gence algorithms like PSO and ABC originally designed 
to solve continuous optimization problems have been suc¬ 
cessfully modified to solve combinatorial problem mm- 
Although SSA only has three parameters besides the pop¬ 
ulation size, it is still very interesting to develop adaptive 
or self-adaptive schemes for SSA to control the parameters 
and reduce the effort in tuning parameters. Last but not 
least, it would be interesting to identify real-world appli¬ 
cations which can be addressed using SSA effectively and 
efficiently. 
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Table 2: Simulation Results for 30-D Problems 





















Table 3: Simulation Results for 10-D Problems 
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